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We present a theory for the spin excitations of the Kondo lattice. We derive an effective Hamil- 
tonian, which describes Fermionic spin 1/2 charge fluctuations interacting with Bosonic triplet spin 
fluctuations. Evaluating the polarization bubble in the magnon Green's function with the free charge 
excitation's Green's function gives a quantitative description of the 'spin gap' and 'charge gap' as 
obtained from numerical calculations for the Kondo insulator. 
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The Kondo lattice model is a generic model of strongly 
correlated electrons. It may be applicable, with some 
variations, to a wide variety of materials, most notably 
the Heavy Fermion metals Jl]]^] and Kondo insulators || . 
It is the purpose of the present manuscript to present a 
conceptually very simple theory for the spin and charge 
excitations of this model, and their mutual interplay. The 
basic idea is, in the spirit of the cell-perturbation theory 
developed by Jefferson and co-workers |J] , to split the sys- 
tems up into sub-units which can be solved exactly. Cou- 
pling the sub-units together the leads to a simple effec- 
tive Hamiltonian. As has been shown in Refs. || for the 
one particle-spectra and as will be shown in the present 
work for the spin excitations, solving this Hamiltonian 
with the simplest approximations possible produces al- 
ready quite satisfactory results. An approach which is 
similar in spirit has been proposed for the Kondo lattice 
by by Perez-Conde et al. and for the description of 
7r-electrons in conjugated polymers by Pleutin |7|. 
The Hamiltonian considered in this work reads 



€k C k<r C k,a + J Y1 S ^,c ■ S n j. 



(1) 



k.cr 



Here £k—J2n ex P{ik-{Rn, — Rm)) t mn is the Fourier trans- 
form of the hopping integral t nm for the c electrons, and 
S n ,c (Snj) denotes the spin operator for conduction elec- 
trons (/-electrons) in cell n. A unit cell contains one o 
orbital and one localized spin. 

The physical parameter regime is J< max \t nm \. It 
was shown previously (H), however, that an 'expansion' 
around the limiting case t mn =0, where the ground state 
for one conduction electron/unit cell is simply a product 
of single-cell singlets, gives in fact a surprisingly good 
description even in the physical parameter regime. The 
basic idea is a mapping between a suitably restricted 
Hilbert space and states of effective Fermions which ba- 
sically stand for charge fluctuations in the singlet back- 
ground. More precisely, when a single cell is in the two- 
electron singlet state, we consider this cell to be empty. 
If the cell n is occupied only by the /-spin we consider 
it occupied by a hole-like Fermion, created by a^ a , if it 
is occupied by three electrons we model this by the pres- 



ence of an electron-like Fermion, created by b\ L(J . So far 
our 'translation table' thus reads: 



\vac) 



a ijvac) -> /,tj0), 



(2) 



In the reduced Hilbert space which is built up from prod- 
uct states of single cell states of the type (||), the Hamil- 
tonian (|l|) takes the form Q 

H ?ff = \ Yl\ ( -£fe + ? T) a k,* a k,<j + Ofe + ^r) & l, A,<x ] 



\ Y si 9n(cr) e k 0L a -/M + H c - 



(3) 



The Fermions have to obey a hard-core constraint, i.e. no 
cell may be occupied by more than one Fermion. This 
constraint is necessary to make sure that the state of any 
given cell be unique. If to simplest approximation we 
relax this constraint, the Hamiltonian (||) can be solved 
by Bogoliubov transformation. Introducing A=3 J/2, the 
quasiparticle dispersion Q is 



£ ± (fe) = (l/2)[e fe ± v / e 2 +A 2 ]) 

with the resulting quasiparticles 

7fe ,-,u = u k b ka + sign{a)v k a^_ ks , 
7h,+,u = -v k b k(y + sign((j)u k a ] _ ks , 



(4) 



(5) 



The choice of signs ensures that e.g. [5 + ,7|]= — 7J.- 
The Hamiltonian (|^) operates only within a reduced 
Hilbert space. Any states where a unit cell with two 
electrons (one / and one c-electron) is in the triplet state 
are projected out. This restriction makes the propaga- 
tion of the charge fluctuations completely coherent, 
and is justified to some degree by the quite good agree- 
ment of the resulting single-particle spectra with Lanc- 
zos diagonalization Q - but ultimately it is of course an 
uncontrolled approximation. Therefore we now want to 
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enlarge the Hilbert space to its full size by introducing an 
additional effective particle £j, whose presence in a given 
cell i implies that this cell is occupied by two electrons in 
the triplet state. Since the three components of a triplet 
are isomorphic to a 3- vector, ti must be a vector particle. 
We thus enlarge the 'translation table' (||) by 



1^0)^2^4(0-^)^/^10), 



(6) 



where er is the vector of Pauli matrices (with cr 2 =3/4) 
and repeated spin indices r, r' are summed over. 
The quantum numbers of the particle are spin 1 and 
charge 0, whence it must be a Boson. As was the case for 
the Fermions, the Bosons have to obey a mutual exclu- 
sion principle, i.e. no two particles, of either Fermionic 
or Bosonic nature, must occupy the same site. 
Having introduced the i ' Bosons to mark spin excited 
cells, we proceed to derive the Hamiltonian which de- 
scribes their interaction with the charged excitations a' 
and W . Figure [l] gives a survey of the different interaction 
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FIG. 1. Different hopping processes which involve triplet 
cells (left panel) and their representation in terms of 
book-keeping Fermions and Bosons (right panel). The tran- 
sition is always accomplished by the hopping of a physical 
f-electron from left to right. Spin and site labels on the 
book-keeping particles in the left panel are suppressed for 
simplicity. 

processes that are possible. The form of the respective 
terms in the Hamiltonian can be inferred from the re- 
quirement of Hermiticity, spin rotation and time rever- 
sal invariance. The first type of interaction process we 
consider involves two Fermion operators and one Boson 
operator. They correspond first to the emission (or ab- 
sorption) of a triplet as a charge fluctuation hops from 
one cell into the other (see Figure |l]a). Another process 
of this type is the pair creation of one electron-like and 



one hole-like charge fluctuation in which a triplet is anni- 
hilated (see Figure |l|b) . The only way to construct a spin 
scalar from two spinors (say: and a) and one vector 
(V) is to first couple the spinors into a vector: 



^mn ~ a \n,T a 'T.T' a n ^ T i , 
^mn b\ nT CTT,T't l n,T'i 



(7) 



and then form its scalar product with the vector t> . Due 
to the product nature of the basis states, the matrix el- 
ements can be evaluated by considering just two cells m 
and n connected by the hopping integral t mn Q. The 
result is: 



Hi — ^ ] t mn [ t m ■ ( S m 

+ (*m — *n) ' ^mn ) ] 



- s: 
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(8) 



Next, we consider processes involving two Fermions and 
two Bosons. Examples are shown in Figure |l|c, where a 
triplet and a charge fluctuation exchange their position, 
and in Figure Qd, where two triplets are annihilated in 
the pair creation of one electron-like and one hole-like 
charge fluctuation. There are two ways to from a spin 
scalar from two spinors and two vectors: first, one can 
contract both the spinors and the vectors separately into 
scalars and multiply these. This term describes processes 
where the Fermion and the Boson exchange their position 
without exchanging spin. Alternatively, one can combine 
the two spinors into a vector following (0) , and form the 
wedge product of the three resulting vectors. We thus 
find the additional terms 
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(9) 



The terms H±, H2 and H3 combined describe all possi- 
ble interactions of the vector Bosons with the Fermionic 
charge fluctuations. The only term still missing to com- 
plete the Hamiltonian is the 'energy of formation' of the 
triplets: H 4 =JY /n t 71 -t n - 

Having written down the Hamiltonian we proceed to in- 
vestigate the physics resulting from it. For simplicity we 
restrict ourselves in this paper to the case of half-filling, 
i.e. the 'Kondo insulator'. As our overall guideline we try 
to keep the approximations involved as simple as possible 
and consequently neglect the 4-particle vertices Hi and 
H3. This leaves us with Hi, which describes the propa- 
gation of a triplet by an RKKY-like process: the triplet 
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may be converted into an electron-hole pair, which may 
later recombine at another site. Alternatively, it may be 
absorbed by an electron-like or hole-like charge fluctua- 
tion, and be re-emitted at a different site. After Fourier 
and Bogoliubov transformation of H±: 
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we introduce the triplet Green's function (where a may 
be x, y or z): 



D(q,u>) = -i(Ttl a t q>a ) 
= lim 



1 



N . (10) 

e^o+ u — J — E(qr, w) + «e 

For the self-energy T,(q,ui) we choose the simplest dia- 
gram possible, i.e. the polarization bubble in Figure 0. 
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FIG. 2. The self-energy £(q,w) used for the calculation of 
the triplet dispersion (the wiggly triplet legs to the right and 
left must be removed). 



(and analogously for the other terms), 
takes the form 



The constraint 
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We now assume that the singlet is condensed Q, so 
that the respective operators can be replaced by the real 
number s which is determined from 



E(<4 



bi 



(12) 



(this means that we neglect the density of triplets, 
tn ' tm which is a good approximation for large and 
moderate values of J) . All in all, in the Hamiltonian (||) 
we then have to replace throughout e k — > s 2 e k , whereas 
H\ changes according to H\ — > sHi . 
Figure || then shows the calculated triplet spectral den- 
sity, A(q,uj)=^(l/n)D(q ) u>). This shows one dominant 
magnetic excitation, wich has quite a substantial disper- 
sion of order J. Such a sharp low-energy mode which 
appears only in the spin correlation function has also 
been observed in Quantum Monte-Carlo simulations at 
low temperatures ||. Near q—n the dispersion takes its 
minimum and thus defines the 'spin gap', A s , i.e. the 
lowest energy required to create a magnetic excitation 
out of the ground state. Also shown in Figure || is the 
dispersion of the two single particle bands, E±(k). Were 
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(11) 



M ' > ~ lo - E_(k) - le + lj - E+(k) + u 

(we are using the free Green's function G M (fe, w') for the 
charge fluctuations 7 T ). 

To describe the low energy scales of the Kondo lattice at 
least approximately, we still have to deal with the con- 
straints. We thereby follow Gopalan et al. |[o| and use 
a simple mean-field like renormalization of the hopping 
integral. We consider the vacuum state, where each site 
is occupied by a singlet, as a condensate of an additional 
Bosonic particle, s^, which stands for to the two-electron 
singlet state. Then, in the Fermionic Hamiltonian (^) and 
in th eterm H\ we have to replace 

Hn,a a n,ij * (^m,<T S «i) ( a n,a S n), 
a ln,a a n,<y * ( a m,<x s m) ( s n a n,a)j 

t ■ S" 1 — > ^ t . <? a 

m ^ ran °n m ^ mn 




FIG. 3. Triplet spectral density A(q, to) for the ID Kondo 
lattice with J/t=1.0, as a function of frequency and momen- 
tum. The Lorentzian broadening e=0.05. 

it not for the low-energy triplet mode, the lowest possi- 
ble excitation would correspond to exciting one Fermion 
from the occupied band into the unoccupied one. The 
energy required for this is the so-called quasiparticle gap 
Agp=i?_(fc=7r)— E + (k=0). If the interaction between 
the particle-hole pair created in this way is small , this 
energy should also correspond to the so-called charge 
gap A c , the energy of the lowest excitation observable in 
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the density correlation function. In fact, density matrix 
renormalization group (DMRG) calculations by Yu and 
White |9) have shown that this neglect of interaction is 
quite a good aproximation. To check our theory, we can 
thus compare the quasiparticle gap Aqp w A c and the 
spin gap, A s , to the results of Yu and White. This is done 
in Figure |I] for different values of the Kondo exchange J. 
The gaps Aqp and A s differ strongly - as is to be ex- 
pected because they correspond to completely different 
types of excitations, i.e. the transfer of a charge fluctua- 
tion from the lower to the upper band on one hand, and 
the propagating triplet on the other. For J/t > 1 the 
agreement between theory and numerics is quantitative 
- this demonstrates that even our very simple and 
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agonalization, as a matter of fact for a wide variety of 
different versions of the Kondo lattice ||. Augmenting 
the formalism by Bosonic spin excitations and using the 
simplest approximations possible, i.e. a mean-field ap- 
proximation to deal with the constraint and a lowest or- 
der evaluation of the Fermionic polarization bubble, we 
could obtain also a quite satisfactory description of the 
spin dynamics. In particular, at least for moderate values 
of the Kondo exchange J, a quantitative description of 
the spin and charge gap as inferred from DMRG calcula- 
tions was possible. This gives reason for some optimism 
to create, mainly by using more sophisticated methods 
to deal with the constraint jfl) in the case of small J/t, 
a theory which is capable of describing the Kondo lattice 
from the largest energy scales (i.e. those of the Hubbard 
bands) down to the smallest ones (i.e. the spin gap). 
It should also be noted that the formalism is indepen- 
dent of dimensionality or lattice geometry. It involves 
comparatively little numerical effort and is therefore not 
restricted either to the somewhat 'toy model type' sys- 
tems considered so far. In fact, both the quasiparticle 
dispersion and the self-energy for the spin excitations in- 
volve only the unhybridized conduction band dispersion 
Using here the results from a 'frozen core' LDA band 
structure calculation may allow to apply the present for- 
malism rather directly to quite realistic systems. 
I would like to thank Dr. C. C. Yu and Dr. S. R. White 
for providing me with the data from Ref. 0. 



FIG. 4. Comparison of spin gap and quasiparticle gap in 
the ID Kondo lattice as obtained by Yu and White M from 
DMRG calculations on a chain with 24 sites to the results of 
the present theory. 

non-selfconsistent calculation of the Green's function ap- 
parently captures the essential physics. For small values 
of J our present approximation is not to be trusted any 
longer. The main reason is that for small J/t the proba- 
bility for charge fluctuations becomes very large (J sets 
the energy scale for the pair-breaking energy) whence our 
mean-field approximation (jl^) predicts a quite dramatic 
renormalization of - such a renormalization, however, 
is not observed in Lanczos calculations ||. Surprisingly 
enough, Figure || shows that despite this problem the es- 
timate for the quasiparticle gap remains accurate even 
for small J/t - this agreement, however, may as well be 
fortuitious. 

In summary, we have derived a theory for the dynamics 
of spin and charge excitations in the Kondo lattice and 
their interaction. The 'purely Fermionic' version, i.e. ex- 
cluding the spin excitations, has recently been shown to 
give quite good results when compared to Lanczos di- 
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